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We study all known and as yet unknown forces between two neutral atoms, modeled as three 
dimensional harmonic oscillators, arising from mutual influences mediated by an electromagnetic 
field but not from their direct interactions. We allow as dynamical variables the center of mass 
motion of the atom, its internal degrees of freedom and the quantum field treated relativistically. 
We adopt the method of nonequilibrium quantum field theory which can provide a first principle, 
systematic and unified description including the intrinsic field fluctuations and induced dipole fluc- 
tuations. The inclusion of self-consistent back-actions makes possible a fully dynamical description 
of these forces valid for general atom motion. In thermal equilibrium we recover the known forces 
- London, van der Waals and Casimir-Polder forces - between neutral atoms in the long-time limit 
but also discover the existence of two new types of interatomic forces. The first, a 'nonequilibrium 
force', arises when the field and atoms are not in thermal equilibrium, and the second, which we 
call an 'entanglement force', originates from the correlations of the internal degrees of freedom of 
entangled atoms. 

PACS numbers: 03.65.Yz, 11.10. Wx, 34.35. +a, 31.15.xk 

INTRODUCTION 

In these notes we delineate a new approach to the study of forces between two (but not limited to two) neutral 
atoms mediated by a quantum (in this case, electromagnetic) field based on nonequilibrium quantum field theory 
PQ. This theory is ostensibly very different from the usual approaches researchers in atomic, molecular and optical 
(AMO) physics are familiar with in the treatment of atomic-optical systems, and it may at first sight appear to be 
too cumbersome or complicated to be necessary. However, with the advances of sophisticated and highly controllable 
experiments in AMO physics made possible by new high-precision instrumentation applied to cold atoms in optical 
lattices (see, e.g., the experiments and the theoretical analysis of [2]) or cavities (with the capability of tracking 
atoms in real time [3]) or nanoelectromechanical systems we are entering an age where traditional theories will 
soon become inadequate. The present method we introduce has the advantage that it is the amalgamation of both 
quantum field theory and nonequilibrium statistical mechanics, the former is required for quantum field (customarily 
referred to as retardation, but there are more involved) effects, the latter for treating processes involving quantum 
dissipation and noises. Not only can this method reproduce all the effects and forces known in the last century as 
detailed below , it can also deal with phenomena and processes more recently brought to central attention from 
quantum foundational and information processing issues, such as quantum decoherence and entanglement dynamics, 
including non-Markovian processes (those carrying memories) which invariably will appear when back-action is taken 
into consideration. Since this method can treat quantum back-action and feedback in a self-consistent manner, it is 
uniquely adept to quantum control considerations [7]. 

Unlike most known treatments of systems in near-equilibrium conditions based on linear or nonlinear response 
theory this is a fully nonequilibrium dynamical description of the atoms' motion [16] . Although our attention is 
focussed on the forces between two neutral atoms derived by finding their relative trajectories (which are determined 
self-consistently by the atoms' interaction with the field and with each other via the field and the field's fluctuations 
in the presence of the atoms), one can also treat radiation reaction, dissipation, and fluctuation phenomena with 
non-Markovian behaviors [5]. As this paper will hopefully illustrate a small initial investment into this new method 
can pay bountifully. 

We consider an assembly of n neutral atoms (labeled by a = 1, ...n) and model the internal degrees of freedom (idf) of 
the ath atom by a three dimensional harmonic oscillator with coordinates Q a , (thus describing the atom's spontaneous 
and stimulated emissions and absorption while interacting with a field). The atoms interact with an electromagnetic 
field (from near-field Coulomb force to far-field radiation) with vector potential A 11 through a dipole interaction, 
but not directly with one another. The force between them arises through field-mediated mutual influences. The 
non-relativistic trajectory of the ath atom is described by z a which, unlike in most previous treatments, is a dynamic 
variable (not prescribed) determined self-consistently by a negotiation amongst all the other variables (Q a ,A^). Our 
interest in this paper is primarily focussed on the center of mass motion of each atom and not on the microscopic 



2 



details of the other variables. The open quantum systems [3] approach can efficiently isolate the desired information 
about the atom's trajectory through a succession of coarse-graining procedures, as detailed below, which take into 
account the overall effects (back-action) of the remaining variables. Using the influence functional method we can 
incorporate the effects of the microscopic physics of the field and the atom's idf and derive an effective equation 
of motion for the atom's trajectory from which the force between the two atoms can be extracted by appealing to 
Newton's Second Law. 

In this setup the dipole moment of the atom modeled by an oscillator is not permanent but only instantaneously 
non-vanishing. The uneven distribution of charge in the atom comes from two effects. First, semi-classically speaking, 
the magnitude and direction of the nuclear-electron separation (which is proportional to the dipole moment of the 
atom) will unpredictably vary in time even in the absence of quantum fields, and second, in the presence of quantum 
fields the atom is polarized by electric field fluctuations. Dipole moment fluctuations are the source of the interaction 
among neutral atoms and are responsible for two types of forces arising from distinctly different physical origin. 

Intrinsic Fluctuation Force 

In the quantum-field conception of a neutral atom the electronic wavefunction surrounding the nucleus has a 
fluctuating component, modeled in our approach by a quantum mechanical harmonic oscillator. As a whole the atom 
will always remain neutral. However, in time intrinsic fluctuations of the oscillator, due to its quantum nature, lead 
to an uneven local distribution of charge in an otherwise (globally) neutral atom which gives rise to an instantaneous 
dipole moment that couples to the attending electromagnetic field. Radiation traveling away from the first (fluctuating) 
atom carries information about the orientation of its dipole (at the time of emission in the past) which eventually 
reaches and polarizes the second atom. The second atom's response to the field leads it also to produce a time- varying 
electric field that travels back to the first atom and is correlated with the activities of the fluctuating atom's idf, 
leading to nonvanishing interaction energy. One can think of the second atom as a transponder which receives a signal 
from the fluctuating atom and then rebroadcasts it. In this analogy the fluctuating atom will receive a signal reflected 
from the transponder atom which encodes its own history. 

This is easily conceptualized if we consider the atoms to be so close that the light transit time between them is 
much greater than all other characteristic time scales governing the dynamics. In such a case the retarded electric 
field is well approximated by the the electrostatic field. Thus, an intrinsic fluctuation of the idf of one atom will source 
a static dipole electric field seen by the second atom. The second atom is polarized by this external field leading it 
too to source a dipole field felt by the fluctuating atom. This process leads to an energetically favorable arrangement 
of the two atom's dipole moments which gives rise to the attractive force between them. 

This type of force due to intrinsic fluctuations in the neutral atoms' dipole moments contain two well-known forces: 
1) the van der Waals force, usually used to describe all interactions between neutral atoms and molecules categorically, 
and 2) the London force which arises from the Coulombic interaction between atoms without permanent multipole 
moments and without the consideration of retardation effects (as Casimir-Polder force does). We refer to forces of 
this type as intrinsic fluctuation forces. 

Induced Dipole Force 

It goes without saying that the quantum field itself possesses intrinsic fluctuations. Any instantaneously generated 
local electric field will induce non-vanishing dipole moments in both atoms. We classify the interaction of dipole 
moments induced by the fluctuations of the quantum field as induced dipole forces. We suggest making a clean 
separation between forces arising from intrinsic (before) and induced (here) fluctuations of the dipole moment of a 
neutral atom because the physical processes produce quite distinct results, as shown in later sections. 

The physical origin of this component of the force is the spatial correlation of field fluctuations. Any given field 
fluctuation will induce correlated dipole moments for the two atoms, much like a long wavelength water wave on the 
ocean will raise and lower two nearby buoys in phase. The excitation of the dipoles by the field will lead to radiation 
that contains information about the emitter. When the radiation from one atom reaches the other the correlation 
between the induced motion of each dipole moment at the time of emission, and subsequent communication of that 
motion via radiation leads to a nonvanishing interaction energy. 

A well known force of this nature is that of Casimir-Polder (CP) [S] who included considerations of the quantum 
nature of the field. This CP force (there is also the CP force between an atom and a mirror which will be treated in 
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our second paper) is a generalization of the London description including retardation corrections as well as effects of 
field quantization - quantization being what imbues the field with its own intrinsic fluctuations. 



For a description of the forces between two atoms we need to know only the averaged effect of the quantum field and 
the oscillator's idf on the atom's trajectory, their details are not of great concern in this quest. Imagine the transition 
amplitude for the total system to evolve from some initial state, \zi n ,ipi n ) to some final state \z out ,(p out ) in time T 
where z labels the atom's position and if is a collective label of the state of all the remaining (environment) variables 
in the total system. Our primary interest is the time development of the atom's center of mass for which the field 
and its interaction with the atom's idf plays a central role through processes like dissipation and radiation reaction. 
For a given final position of the center of mass there can be many consistent final field and oscillator states, likely 
unobservable. Summing over all final environment states compatible with the atom's motion is necessary when we are 
ignorant of the final state of the environment whether we choose to ignore those details or they are not measurable. 
This leads to an effective transition amplitude for the trajectory of the atom alone where all environmental effects 
on the trajectory have been taken into account. Carrying out this process of coarse — graining where the final field 
and oscillator states are traced over leads to an effective action that self-consistently accounts for all back-action of 
the field and the atom's idf on the atom's trajectory. The equation of motion for the atom, and thus the atom-atom 
force can be obtained through a variation of this action [8]. 

As we shall show the present formulation goes beyond previous work in that we can derive the forces between two 
atoms for fully dynamical and under nonequilibrium conditions. When the spacing between the atoms is held fixed we 
recover the well-known London and CP forces. For the case where the atoms and field are not in thermal equilibrium 
we find a novel far field scaling for the induced dipole force diminishing as l/z 3 rather than l/z & , and for the case 
when the two atom's are entangled we find a novel near field scaling that enters at second order in perturbation theory 
as q 2 /z 2 as opposed to the standard q^/z 7 where z quantifies the interatomic distance and q the electronic charge. 
To the best of our knowledge this nonequilibrium force and the entanglement force behavior have not been reported 
in the literature. 

This paper is organized as follows: In Sec. 2 we introduce the microscopic details of the system by defining the 
action describing the dynamics of the entire system. The worldline influence functional method is adopted where 
the environment degrees of freedom (field + oscillator) are traced over to find the time evolution of the reduced 
density matrix. The equations of motion for the atomic trajectories are then obtained from the saddle points of the 
reduced density matrix. In Sec. 3 the explicit form for the atom-atom force is computed, and physical origin of each 
component is explained. In Sec. 4 an initially entangled state for the two oscillators is considered. It is found that 
the correlation among the oscillator coordinates leads to a contribution to the atom-atom force that enters at second 
order in perturbation theory. Sec. 5 concerns the possibility of detecting the nonequilibrium atom-atom force (when 
the field and atoms are not in thermal equilibrium) which is new and the novel entanglement force we discovered as 
outlined in Sec. 4. 



Coarse-graining and Back-action 



THE MODEL 



We describe the microphysical degrees of freedom of the entire system through the following action 



S[Q a ,z a ,A»] = Y,( s QiQa] + Sz[z a ] + S mt [Q ai z ai A»]) + S E {A»] 



(1) 



a 



where the sum is over all atoms. The action describing the dynamics of the oscillator is given by 




(2) 



where [i a is the ath oscillator's reduced mass and A its worldline parameter, Q a being its natural frequency. The 
electomagnetic field action is given by 
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S E [A»]=- l -j ' d^xF^F^ (3) 

(the subscript E stands for the electric field) where is the 4-vector potential and F^ — d^A u — d^A^ is the field 
strength tensor. The action for the center of mass motion is 

(4) 

where M a is the atom's total mass and V[z a ] is an external potential. 

In the dipole approximation, the potential energy for an atom interacting with the photon field takes the form 
— q Q ■ E[z\, where qQ is the atom's instantaneous dipole moment and E is the electric field leading to the interaction 
action Si n t[Qa, z a , A^] = q a J d\Q l a (\)Ei[z%(\)]. Above, q a , quantifies the coupling of the ath atom to the field. 
[Greek indices will refer to spacetime components of a four-vector, zero referring to time, and Roman indices refer 
to spatial components where we will exclusively use the letters {i,j, k} to avoid confusion with the letter a used to 
label atoms. Contraction of four-vectors is undertaken with the Minkowski metric with (-,+,+,+) signature, and the 
Einstein summation convention is used throughout. 



Sz\z a ] 



d\ ^z a 2 (A) - V[z a ] 



Worldhne Influence Functional 



Assume that at time ij n the quantum statistical state of the oscillators, trajectory and field is described by a density 
operator p(ti n ). This state is unitarily evolved from the initial time ti n to a later time t > ti n , and can be expressed 
in terms of path integrals by considering matrix elements in an appropriate basis. 

To isolate the influence of the field on the dynamics of the atom we coarse-grain over the field variables to construct 
the field-reduced density matrix, p r (Q a ,Q' a ; z a , ^t) = J dA^ p(Q a ,Q' a ; z a , z! a ; A 11 , A^;t). By assuming that the field 
is initially uncorrelated with the other degrees of freedom the reduced density matrix takes the form, 



Pr (Q a , $ a ;z a ,z' a ;t) = YlJ dQ m , a dQ' m a J dz in , a d? m a VQ a VQ' a jf " Vz a jf ' V* a 

X g ( q[ a] z[ a] q[ al Z ^' a ^PQa(^Qin,a^Qi n .aT^in)Pz{zi nai Z inal ti n )F\ L J^ ,J ] (5) 

which introduces the influence functional (IF) F\J^~ , J v+ ] [TT]. If the initial state of the field is Gaussian in field 
variables (which includes vacuum and thermal states) the influence functional can be calculated exactly for the dipole 
field interaction. 

F[Ji*-,r+] = expjz J d*y J d i y'[J^(y)D^(y,y')r+( y ') + l -J»~ (y)D^(y, y')J"~ (y>)} } (6) 
Here the current density is 



J u (x) = -$>«/ d\ Kifi 5 A (x» ~ <(A))<X(A), (7) 

J + = ( J + J')/2 and J~ = J — J' are its semi-sum and difference, respectively, where prime distinguishes histories, 
and Ki^ = dirjQft — d^rji^ is a differential operator that relates the photon field to the electric field by contraction 
i.e. E % = k^A 1 *. D r ^ e J(y,y') and D^{y,y') are the retarded Green's function and Hadamard function for the field 
respectively. In the Feynman gauge they can be expressed in terms of the retarded, D ret , and Hadamard, Dh, Green's 
function for a massles scalar field. 

D r £{x,x') = V ^D ret (x,x') D^(x,x') = ^Dnix.x') (8) 

At zero temperature they take on the explicit form 

D ret (x,x') = ~6(t-t')6(a) D H (x,x') = --^ r (9) 

where a is Synge's worldfunction defined to be half the geodesic distance between the four- vectors x and x' , a = 
{x-x'f/2. 
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Oscillator-Reduced Influence Functional 

We isolate the net influence that the oscillator's idf Q a and the field have on the trajectory by tracing over 
all final oscillator configurations p or {z a ,^ a ]t) = J dQ a p r {Q a ,Q a ;z a ,z' a ;t) introducing the oscillator-reduced density 
matrix p or . 



(10) 



All the effects of the environment are now packaged in the oscillator-reduced IF, Fz\z~ , z + ). The development has 
been simplified by working in the rest frame of the second atom in so doing z-i is no longer treated as a dynamical 
variable. 



F Z [H,4]=T\ [ dQ a dQ mM dQ' in . a VQ a C^e^^^l-^^^pQ^ftn.^Q-n.a;^)^ 7 ^",^] (H) 

a J JQin.a J Q'i„,a 



To elucidate our approach we write ( 11 ) in a more suggestive form 



5ji ' ^ Jo 



(12) 



j.±=o 



which defines the IF for a three dimensional harmonic oscillator, F [j+,j a ].To bring F[J^ , J v+ ] out of the 
path integrals in (11) [Q„ (A)]™ is replaced with functional derivatives on the IF for the harmonic oscillators 

. The explicit form for F a [j + , j~ } is 



j»±=o 



F a {3 + J-} = I dQ a J dQ in , a dQ in , a p Qa (Qin,a,Qin, a ;tin) L' VQ a j / VQ'J 3 * tf«M S «J 

(13) 

In the above the dot product between current and oscillator coordinate i.e. j ■ Q is taken with respect to a three 



dimensional Euclidean metric. For a Gaussian initial state ( 13 ) can be evaluated exactly 



FaUiJa] = -^exp i d\dX'[j-(X) ■ j + (AV et ,a(A, A') + -j~(X) ■ j~ (X')g H>a (X, A')] 



(14) 



where g r et,a{X, A') and gu,a(X, A') (expressed below at T = 0) are the retarded and Hadamard Green's functions for 
a one dimensional harmonic oscillator with natural frequency il a , mass p a , and N is a normalization constant. 



9ret,a(X, X ) 



0(A-A')sinfi o (A-A') g B , a (\,\') 



1 



COS Sl a (A — A') 



(15) 



Decoherence and the Semi-Classical Limit 



The complex norm of ( 11 ) at leading order in a ^"-expansion follows 



\F Z \Z[,4]\ = cxp 1^- J dXdX' z\- 



(X)N ZJ (X,X')z{-(X>) 



(16) 



where is a symmetric positive definite kernel. Thus, we observe that the off-diagonal elements of the density 
matrix in ( 5| are strongly suppressed for large values of z~ = z — z 1 as is indicative of decoherence of the quantum 
trajectory |12j . Decoherence of the trajectory due to its interactions with the quantum fluctuations of the environment 



G 



and the internal degrees of freedom of the atoms permits the existence of a semi-classical limit for the oscillator's path 



through space. Using a saddle-point approximation to evaluate ( 10 ) one can show that the semi-classical dynamics is 
determined from the variation 



SS C gea[z? \z t ] 



6zt(r) 



= 0=^ Mz k {j) = A(r) (17) 



where the so-called coarse-grained effective action is given by Scgea[z\ , z\~] = 5z[£i] — S*z[^] + Sif[z^ , z[], and 
Sif[zi,Zi] — —i^nJ~z[z[,z^] defines the influence action. The force acting on the trajectory due to its interactions 
with the oscillators and field is given by 



5Si F [z^,z 1 

Jk\T) - 



6z?-(t) 



(18) 



For general atom motion this force contains all known effects, including the Lamb shift, radiation reaction, dissipation, 
and the atom-atom force. 

NONEQUILIBRIUM ATOM-ATOM FORCE 



The suppression of the reduced density matrix for off-diagonal elements justifies an expansion of ( 12 ) for small 
values of z*[ . The linear order term yields the influence force and is represented by an infinite series in powers of the 
coupling. The local (spatially independent) terms in this expansion lead to the aforementioned Lamb shift, radiation 
reaction, and dissipation. The atom-atom force can be obtained from this series by extracting the terms that depend 
upon the spatial separation of the atoms. 

To simplify the presentation we rewrite the influence functional for the atom's trajectory as Tz = (e lScff ) where 
the form of S e ff = — iln T can be taken from ^ 

S eff = J d*xj d i x'[J^-(x)D; e J(x,x')r + (x') + l l J^-{x)D^ v {x,x , )J v -{x')] (19) 

and (...) is the expectation value with respect to both oscillators under the condition of no interactions. 



(->o = n/ ' d ^ d Qin,ad& in , a PQAQin, a ,QL,a,Un) J^" VQ a PQ^^-So [<3' a D (...) ( 20 ) 



S e f f is a quadratic function of the current density ^ which depends on a sum of delta functions with support at 
each atom's position. Thus, one can see that the cross terms S cross between the two atom's currents appearing in 
S e ff will lead to terms that depend upon their spatial separation. 

S cross = [ d^xyrWD^^x'^+ix') + J^-{x)D^{x,x')J^{x') + ljr(x)D"(x,x')J^-(x')} (21) 



J£ refers to the current density of the ath atom, with a — 1 referring to the distinguished atom, where all the 
atom-atom forces we are studying here act upon. {S cross ) , the expectation value of S cross , will vanish for initially 
uncorrelated and Gaussian oscillator states because it is linear in the coordinate of each oscillator. Therefore the 
leading order contribution to the atom-atom force will be proportional to the square of S crosa at order qfq%. Expanding 
Tz in powers of S e f / we express the IF as 

T z - e iSl ^ w 1 + (self energy terms) - * ((S cross [zf, Q± ]) 2 ) q . (22) 

The leading order linear terms O(q^) contain the back-action of the field on the motion of the atom itself only. We 
refer to them as "self energy terms", borrowing a terminology from particle physics, and for a stationary atom these 
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effects are unimportant. We focus on the quadratic term which contains not only higher order self energy type effects 
but also the leading order contribution to the atom-atom force contained in S cross . 

Note that S e ft contains a term linear in the retarded propagator for the field. It makes sense that the force will 
manifest as the square of (21) because, as described in the introduction, fluctuations in the dipole moment of one 
atom induce radiation that travels to the other, influences its dynamics, and induces the other atom to radiate. From 
a diagrammatic viewpoint this process requires two propagators. 

The leading order expression for the force can be derived from 



fk(r) 



SSif [zf] 



5z 



k — 



2 5z k r(T) 



(23) 



where the first equality holds for a stationary trajectory (radiation reaction and dissipation vanish). 
Carrying out the variation in ( 23 1 yields three contributions the first two being 



fk( T ) = 



dX 



Ai 



d A x 



d A yg r etAr, \)E r { t {z a (\),y)G H {x,y)d k {x')Elf{x',x)\ x , =za{T) 



(24) 



fk ( T ) = 2«i«2 



dX 



d A x I d 4 y gnAi 



X)E r { t (z a (X),y)G ret (x,y)d k (x')E^(x\x)\ x , 



«(r) 



(25) 



where E i j(x 1 x') is the dyadic electric field Green's function E i j(x 1 x') = Tr{pEEi(x)Ej(x')} which arises in our 
formalism through a contraction of the operator k 1 ^ with the photon field. Note also that after all functional derivatives 
are taken z\ — > z, G re t{x,x') and Gh{x,x') are meant to generally represent the retarded and Hadamard functions 
for whatever the atom interacts with. For example to find the surface-atom force the integration volume V is taken 
to be the half space, and the Green's functions describing the physics of the media occupying that region are used. 
For the specific case of an atom located at z 2 , G(x,x') cx g 2 (t, t')5 3 (x — z 2 )(5 3 (x' — z 2 ). More general cases will be 
considered in a future paper. 



The form of (24) and (25) can be explained by appealing to the heuristic description of the force given in the 
introduction. f A and f B arise from the intrinsic fluctuations in the dipole moments of the atoms. This can be seen 
by noting that they contain the atom's Hadamard function i.e. the symmetric two point function for the oscillator 
degree of freedom. The two retarded electric field Green's functions account for the transfer of information between 
the two atoms, and the retarded Green's function for the atom characterizes its response to an external field see Fig[l]. 

The third component of the force, / arises from induced fluctuations of the atom's dipole moments. The 
retarded Green's functions for the two oscillators, g re t,a, characterize their response to a given field fluctu- 
ation. The kth component of the induced dipole moment of the ath atom can be written as d k nd = 
q a J dX' g r et,a{\ ^')E k [z"(X')] where E k [z"(X')] is the fcth component of the electric field at the position of the atom. 



The symmetric two-point function of the induced dipole moment quantifies its fluctuations \{d\ nd a (t), d\ nd b (t')} 

QaQb J dXdX' g re t,a(t, X)g rett b{t' , X')E^ H (z" (A), z"(A)). The remaining electric field propagator, EJ^, carries information 
about the motion of one atom to the other and accounts for the form of f c see Fig[l] 



Kir) = 2<fll2 



dX I d 4 x I d 4 y g ret AT,X)[d k (x')E r l t (x\x)G ret (x 1 y)E^z a (X),y)\ x , =za{T) 

+E% t (z a (t) , x) G ret (x, y)8 k (x')E^ (x' , y) | X - = ^ (A )] 



(26) 



where d k (x) denotes differentiation with respect to x k . The previous form of the force is valid for any atomic motion. 
However, a self consistent treatment would require that the aforementioned 'self energy' terms be included in order 
to account for the back-action of the field on the atom itself. 




FIG. 1: The illustrations depict the physical origin of the intrinsic fluctuation and induced dipole forces. On the left intrinsic 
dipole fluctuations (represented by shaded oval); 1. radiate information about their motion, and 2. this radiation induces a 
correlated dipole moment in the second atom (solid black arrow denotes an induced dipole moment). The induced motion at ti 
leads to radiation that travels back to the fluctuating atom. At ts the radiation produced at step 2 will produce a local electric 
field near the fluctuating atom which carries information about its own fluctuations in the past. The illustration on the right 
depicts the physical origin of the second component of the force arising from field fluctuations and their spatial correlation. Step 
1 shows how a field fluctuation induces correlated dipole moments in both atoms. The induced motion of the dipole moments 
will lead to radiation emitted from both atoms containing information about their motion (only left moving radiation included) . 
At t2 the radiation generated by the induced motion produces a local electric field around each atom that is correlated with 
its motion. 



Induced Dipole Force 



In this section we calculate the induced dipole force explicitly by plugging in the retarded Green's function for the 
second oscillator, g re t,2, and choosing z 2 to be the origin of coordinates. 

fk{r)= l -qlq 2 2 ' dX J dt J dt 1 g retA ( T , A) flpet)2 (t, t')[E% t (z a (r), t' , d)d k (x')E?](x' ,t,d)\x>=z° W 

+E%(z a (X),t', 6)d k (x')Elf(x', t, 0)\ x , =za{T) ] (27) 

The derivatives operating on the various Green's functions can be simplified by employing the equation of motion for 
the field, the resultant form valid for general atom motion follows. 
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fk ' = qj ^ dX dt dt'g retA (T,\)g reti2 (t,ty(T 1 t')\s"'[a(z a (T);t' 1 0)} 



-D H [a(z a (X);t,6)} 



xa(z a ( T ); t', 0) k ( <r(z a (X); t, 0), a(z a (X);t, 0) J a(z a (T);t', 0), a(z a ( T );t', 0) 4 + [a(z a (X); t, 0), a(z a ( T ); t', 0))f 

'1 



-S [a(z a (r);t',0)] 



-D H [a(z a (X);t,0)} 



xa(z a (X);t,0) k \a(z a (X);t,0) j a(z a (X);t,6y <T(z a (r);t', 0); a(z a (r); t', O) 4 + [a(z a (X); t, 6), a(z a (r); i', 0))f 
+4W« a (A); t, 0)] (&a{z a {T); t', 0) t a(z a (r); t', 0)V(z tt (A); i, 0) fe + 2 ( r(z Q (A); t, 0), a(z Q (r); t', 0))' a(z a (r);t', 0) fe 



+<5V(z Q (r);t',0)] 



4 Djj [cr(z a (X); t, 0)]a(z a (T);t' , 0)j a{z a (r);t' ,6)* a(z Q (A); 0) fc 



+20 £> ff [a(z a (A); i, 0)]a(z Q (A); t, 0) fe 



(28) 



Here, primes on functions denote derivatives with respect to a, and Ofc = d k u denotes differentiation of a with respect 
to x k . 



We can separate ( 28 ) into 4 terms with differing number of cr-derivatives and specify a static trajectory for the 
distinguished atom to bring the derivatives outside of the integral i.e. d/da = z~ x djdz. To distinguish which Green's 
function a given cr-derivative acts on we attach a dummy subscript to z that should not be confused with an atom 
label. Once all derivatives are taken Z\ and z 2 are set to z, the separation between the two atoms. The evaluation of 
the t-integral can be done by substituting 8{a{x,x')) — — ~^~^ X ^ ■ 



jjfgg 1 
4tt Z 2 



1 d \ 3 / 1 d x " 



Z2 dz2 J \ Z\ dz 



1 d \ 2 ( 1 d x ; 



Z2 dZ2 



zi dzi 



I 



dX / dt g re t,i( T > ^)9ret,2(t,T - Zi)D H [cr(z% (X),t,6)] 



2z q 



(29) 
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mi 

47T 



2 1 

mi 

47T 



1 d \ 2 f I d V 1 



Z2 dz2 



z\ dz\ 



Zl 





-V 




d \ 






dz 2 J 




dz\ j 





dX j dt g re t.i(T, X)g rett i{t,T - Zi)D H [a(z% (X),t, 0)] 
dX / dt g ret ,i(T, X)g ret .2{t, r- zi)D H [a{z% (X),t, 6)] 



8z 2 



4z' 



(30) 



(31) 
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47T 





-V 



















— dX dt g ret ,i(T, X)g reta {t,T - z-i)D H [<j{z% (A),t,0)] 
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(32) 



We can express the Green's function for the field through a mode sum and subsequently evaluate the A and t 
integrals in the long time limit. The exact field-influenced dynamics of the oscillators will be dissipative, however this 
dissipative effect does not appear at this order in perturbation theory, but can be modeled phenomenologically by 
inclusion of an infinitesimal dissipation in the oscillator equation of motion i.e. g re t,a{t ~ O — ^ 9ret,a{t ~ £')e _e ('~' '. 

At finite temperature the Hadamard function for the field can be obtained through periodicity in imaginary time, 
or by taking the trace of symmetrized field operators with respect to a thermal density matrix. Expressing the result 
as a mode sum we find. 



f-Cl 
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2 1 

- — z - 

7T 2 



1 d 

Z2 dz 2 



1 d 

Zi dz\ 



\ d 

Z2 dz 2 



I d 

Zi dz\ 



Z\Z2 Jo 



duj ai(oj)a.2(u)) coth(/3o;/2) $m.wz2 cosw^i 



2z q 



(33) 
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fC2 



2 / 1 d 

7T \Z 2 dz 2 



1 d 

z\ dz\ 



ziz 2 Jq 



dui ai(uj)a2(oj) coth(/3w/2) sinwz2 coswzi 



(34) 



fC3 



1 d 

Z 2 dz2 



1_ _d_ 
Z\ dz\ 



Z\Z2 Jq 



du a.\(uj)oi2{<jj) coth(/3w/2) sinwz2 cos uizi 



4z< 



(35) 
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1 d 

z 2 dz 2 



1 d 

z\ dz\ 



1 



Z\Z2 Jo 



dui ai(ui)a 2 (oj) coth(/3w/2) sinwz2 coswzi 
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(36) 



As the retarded Green's functions for the atoms characterizes the response of their dipole moments to an external 
field they play the role of the dynamic polarizability, a. The result above is expressed in terms of the frequency 
dependent form, a a (tu) — g^(47r/i o (0^ — w 2 )) -1 , which can be derived from the classical equations of motion (the 
aforementioned infinitesimal dissipative term kills an imaginary part in the infinite time limit). (3 is the field's inverse 
temperature. All the derivatives can be taken and the force can be expressed as an integral over frequency. 



/f = 



duj ai(w)a2(w) coth(/3w/2) 



wz(18 - 8z ui + z u) ) cos 2ujz 



+ (-9 + 16wV - 3wV) sin 2wz 



This expression agrees with what can be found in the literature for the CP force in a finite temperature field 
the far-field at zero temperature we recover the well known form [14] 



(37) 



In 



f° «-^ax(0)a 2 (0)l 



(38) 



where a UV regulator must be employed to render the frequency integrals finite. 

If however we take the dissipation to be zero (as it truly is in our perturbative approach) the force is altered because 
the polarizability acquires an imaginary part a(uj) —> (g 2 /47r/iS7)[17/(17 2 — w 2 ) -1 + i(ir/2)5(uj — 17) — i(ir/2)5(uj + 17)]. 
The imaginary term plays an important role when the quantum nature of the dipole moment of the oscillator is 
accounted for. When such a term is neglected, the contribution to the atom-atom force from f A and f B dominates 
in the far field as 1/z 3 rather than 1/z 8 at T = 0. We denote this contribution to the force by 5fg. 



5f z 



2 2 



167T 2 /ii^ 2 17il7 2 z 7 Jo 



dco coth(/3w/2) 



17i5(cj - 17 2 ) 



+ (17i <-> 17 2 ) 



9 - 2z z lo z - z*lj* + (9 - 16z z uj z + 3z 4 w 4 ) cos(2u;z) + zu(18 - 8z z uj 2 + z 4 cj 4 ) sin(2u;z) 



As z — > oo this term possesses the asymptotic scaling, 



(39) 



2 2 



167T 2 ^i^ 2 ^1^2^ 3 Jo 



dcu lu 4 coth(/3w/2) 



17i%-Q 2 ) 



(fii o 17 2 ) 



(40) 



but at T = for z —> 6f is subleading to the dominant 1/z 7 near field scaling from the London term. 



Intrinsic Fluctuation Force 



The treatment by London of the atom-atom force can be reproduced by computing the interaction energy of two 
atoms interacting via the Coulomb potential. The force follows from the negative gradient of the perturbed energy 
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eigenvalues. We obtain an analogous expression for the London force in our formulation but with an additional 
contribution from retardation effects as we treat the held relativistically. 

The contributions to the force from and ff can be computed in the same way as the contribution , and so 
we omit the details of that calculation here, and only state the result in the long-time zero temperature limit. 



167T 2 /ii/i2^1^2 ^1 



01' 



2V&z 2 



(41) 



fz can be obtained from ff by exchanging fii and f^- These terms are responsible for the near-field behavior and 
agree with those derived by London when retardation corrections to the field Green's function are neglected [4]. 



Lon 



9<Zi<?2 



167T 2 fiifj,2^ll^l2 ^1 + ^2 Z 1 



(42) 



The thermal version of the previous result does not make sense for a single oscillator where temperature is an 
ill-defined quantity, but does in the case of a gas of atoms. If the gas is sufficiently dilute the force between two 
collections of trapped atoms can be approximated using the density distribution of the gas and f z |15j . The finite 
temperature form follows where (3 a is the inverse temperature of the ath oscillator (or trapped gas). 



Si 



1 2 

qfqi 



{l[ [9 + 2Qjz 2 + fljz 4 } coth(#2fi2/2)^7 



167T 2 /I i /i 2 1^2 ^1 — ^2 

In the far field the leading order behavior reduces to the following form 



mi 



coth(/? 2 2 /2)^ 



(43) 



(44) 



Note that when the field and the atoms arc in thermal equilibrium this new asymptotic scaling cancels with an equal 
and opposite contribution contained in Sf^ and the standard far field scaling 1/z 8 is restored. When the atoms and 
field are out of thermal equilibrium this cancelation no longer occurs and the dominant contribution to the force scales 
like 1/z 3 , where the zero temperature contribution cancels as indicated below: 



2 2 



87T 2 jLti/i2 £l\ 



01 



(45) 



ENTANGLEMENT FORCE 

The previous derivation of the atom-atom force assumes that the initial state of the two oscillators is uncorrelated. 
If however, the two atoms are initially entangled then a new contribution to atom-atom force arises. To our knowledge 
this force has not been reported in the literature. 

We begin by computing the oscillator-reduced IF for two initially entangled atoms. 



^z[*r,#] = n 



Writing ( 46 ) in a more suggestive form 



dQ a dQ in ,adQ' ma VQ a pg;^ [<3'J) 

XpQ(Qin,l, Qin,li Qin,2, Qin,2i *in)<7V , J 



(46) 



z a ,z a , * * 

°Ja UJa 



F[it iji lit 1 32 



(47) 
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defines the IF for two entangled harmonic oscillators, F[j^~ , j 1 , j 2 , j 2 }■ To bring T[J^ , J v+ ] out of the path integrals 
in ( 46 ) we replace the oscillator coordinates with functional derivatives as before. 



1 = II y d Qa J AQin.aAQ in, a PQ^Qin.li Qin,li Qin,2i Qin,2i 



Qa , _ _ _ 



For the initially entangled squeezed Gaussian state 



(48) 



PQiQin.li Qin,l> Qin,2i Qin,2> ^in) 



it a 



exp 



P 2 ((Q in ,i + Q m ,2? + (Qi„,i + QL2) 2 

+ ((Qin,l - Qin,2) 2 + (Qm,l ~ Q'in.2) 2 

or \ 



(49) 



(48) can be evaluated exactly. For this case, like oscillator coordinate components are entangled together with equal 



magnitude in each direction i.e. the parameters a and j3 are common to each component. 
The influence functional for two entangled oscillators follows 



J2J2) = -Afexp I i UB A flfA + *B s gs - o?B\ - - ^g\ - /3 2 .g| + iip 

where Af is a normalization constant and the definitions of i? Q , g a and ip are defined below. 



(50) 
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B a = J dt J~(t)[cotD,Tsmn(t-ti)+cscQTsmQ(t f -t)] 
g a = ^ J dt J~(t) sinfi(t - ti) 
E J dt 9aJ a (t)smn(t f -t)-^n g 2 a cotnT+ l - J dtdt'{J a (t)J a (t')gUt,t') + j' a (t)j:(t')gUt,t')} 



(51) 



(52) 



(53) 



The subscript £ denotes Cs = Ci + C 2 , and the subscript A denotes Ca = C\ — C 2 ■ 

The entanglement force comes from the leading order contribution to Tz that depends upon the spatial separation 
between the atoms. Previously we needed to consider the square of S cross . However, when the two atoms are entangled 
there exists nonvanishing cross correlation between their coordinates such that (S cross ) =/= 0. 

So in distinction to the previous section we have 



jr z = e ^i^j ~ ! + (sclf energy terms) + z ( 5 cross [^, Q±] ) + 0(0 
where the force can be derived from 



(54) 



fk nt (r) 



1 _ . \ S cross [z a , Q a 1 



(55) 



Expanding S cross for small 2* we arrive at 



5 + g l92 / d\d\' Zi~ (X)[dkKi Kj D^(z? + (A), z 2 a (A'))Ql + (A)Q J 2 +(A') 



+jd k >4 / t p- t (^(A),^ + (A'))Qr(A)Qr(A') + ^ fe <<'^(^ + (A),^(A'))Ql + (A)Qr(A')] + 0(O (56) 
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where a prime in the index of a derivative operator means differentiation with respect to the second argument. Only 
one term survives after we take the expectation value, that which contains the cross correlator (q 1 + (\)Q j 2 + (A')) 
which equals 



(Q\ + {t)Qi + {t')) o = & g e nt(t,t>) = i -/3 2 ) & 



1 



■ sinf2(t — ti) sinf2(t' — i,) - 



esc 2 nTS(t)S(t') 



(57) 



(fin) 2 v lJ y 11 4/3 2 

where T = tf —t{ and S(t) = sin Vt(t — tf) — sinQ(t — ti + T). After taking the expectation value and then using (55 ) 
we obtain the entanglement force. 



fk n \r) = - gi q 2 J dX g ent (r, X)d k ^ K iv D r * (*? (r) , z«(A)) (58) 

All derivatives can be taken on the field's retarded Green's function and simplified using the equation of motion. We 
then specify the trajectory to be static to arrive at 



5z 



Dret(cr). 



(59) 



With a static trajectory specified the a derivatives can be expressed in terms of z derivatives and can be factored 
out of the integral. The retarded Green's function can be expressed as a delta function D ret (cr(zi(r), z 2 (A))) = 
5(t — A — z)/(A'Kz). We work in a coordinate system centered at the a = 2 atom, where various forces act upon (also 
the origin of the xy-plane), with z axis along the ray connecting the two atoms at distance z apart (pointing from 
atom 2 to atom 1). Thus the distinguished atom (a = 1) is located at (0,0, z). This leads to the explicit expression 
for the entanglement force. 



f! nt (r) = 



<M2 

47T 



1 d 

z dz 



In the infinite time limit r — > ti 



if ->■ oo g ent (T,T- z) - |(/3 2 - l/a 2 )(a 2 //3 
in the far field but has a well-defined near field limit i.e. £lz 0. 



+ 5z 



ld_ 

z dz 



1 



z 

2 lal 



{ (r,r- z) 



(60) 



l/(fitt) 2 ) cos Slz. The force vanishes 



n 



<M2 

32vr 



1 



1 



,2 



(61) 



J 2 (/xO)^ 

This effect is not only due to entanglement between the two atoms, but it is also due to retardation. For the case 
considered above where the degree of entanglement between like components of the atom's dipole moments has the 
same magnitude the interaction energy as described through the Coulomb potential vanishes. Thus, only through the 
inclusion of relativistic effects does any force manifest. 

However, the previous discussion can easily be generalized to the case where the magnitude of the parameters a 
and j3 is not common to all directions. As such the initial state for the oscillators takes the generalized form 



PQ 



{Qin.l, Qin.lt Qin.1t Qin,2t tin) — J][ 
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\irajj { 



% ((Qi.x + QLaf + (QL,i + QL,2f 



+^((Qki-Q 



J \2 
in, 21 



{Qin.l 



(62) 



The development follows closely that given for the previous case and so we only state the result in the near-field 
long-time limit, 



in 



9l<?2 



A, 



A y - 2A 2 



1 



(63) 



where we have used the shorthand A,- = |(/?| — l/ a |)( Q; f /Pj ~ Vl/ 1 ^) 2 )- Note that if the parameters a and /? are 



equal for all directions (63) vanishes. The sign of the force can also be changed by the appropriate choice of the 



squeeze parameters a and j3. 
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POSSIBILITY OF DETECTION 
Atom and Field Out of Thermal Equilibrium 

In this section we compute the relative magnitude for the atom-atom force when the field and atoms are not in 
thermal equilibrium to the force at zero temperature. We focus our attention on the case where the atom's are in 
their ground state and the field is in a thermal state of inverse temperature ft. Measuring this new asymptotic scaling 
requires a balance between temperature and the first optical resonance of the atomic species used. For the case when 



£1/3 >> 1 ( 45 ) is exponentially suppressed, this would rule out the use of heavier atoms like Rb near room temperature, 
the only hope is to work in the regime where f3£l > 1, not only to prevent suppression by the Planck factor but also 
to prevent the excitation of the atom so that measurement can be done before thermalization. 



The relative magnitude of (45) to fc in the far field shows when this new scaling will dominate. For realistic 
experiments, atom-atom distance of the order of /im, the high temperature limit is beyond access for the temperatures 
and atomic species we are considering, so we replace 6f c with its zero temperature form (also the appropriate factors 



of c have been restored to ensure that ( 64 ) is dimensionless) 



8f c 8tt n?^5 1 

^ * _*o, "a) (64) 



f c i6ic 5 nf - nl - 1 

If the atomic species are the same the previous expression reduces to 



f c 



24tt fQz 
322 



1 



Tuning to hydrogen's first optical resonance (fi ~ lOeV, SI « 2.4 x 10 15 Hz, or fl ps 116, 000K ) we find 

5 



f c 



24tt / 8 



322 \fim ) e 



1 



m _ i ' 



(65) 



(66) 



If the atomic species are different we find different behavior. Particularly, when one of the atom's first optical 
resonance is very large such that /3 ft >> 1 (like Rb near room temperature) the Planck factor for that atom will be 
strongly suppressed so its contribution to the force can be ignored. In such a case (|64|) takes the form 



6f c _ 8^ fftz\ 5 8tt (8z\ 5 1 

f c 161 ^ c J e^-l m\iim) e^-l [ ' 

where we have a different sign and a slightly different coefficient. 

For atoms Sf c / f c only becomes significantly greater than 1 for large distances and very high temperatures and so 
is unlikely observable. However, these effects may play a role in the laboratory for molecules with sub eV excitation 
energies. We leave a study of those effects for a later work. 



Entanglement Force 



Now that we have an expression for the entanglement force at short distances we check for regimes in which ( 55 ) 
will dominate. To do this we take the ratio of the entanglement force to the near-field van der Waals force. After 
restoring all physical constants to yield the correct dimensions and allowing both atoms to be the same species we 
find. 



jLon g c 2 q 2 \r ^ 2 J V/3 2 

Above /3 — /3/fiQ, and a — fifla. 



^^-i)f|-.|^ (68) 
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FIG. 2: Relative magnitude of the entanglement force to the London force. Intersection of the various graphs with the black 
line yields the value of the interato mic d istance where the the two forces are equal in magnitude. The various colors represent 
different values for the prefactor in (69 1 i.e. P — 8.9 x (ji 2 — — 1^. For red P — 0.01 for blue P = 0.1 and for green 

P = 1. The plot shows that for interatomic distances satisfying the near field condition flz « 1 that the entanglement force 
dominates for order 1 prefactor for distances of a few nanometers. 



By tuning the frequency to the first optical resonance of Hydrogen, taking the reduced mass to be the electron mass 
and q to be the electronic charge we find. 



f L 



8.9 x LS 



ft 2 



(-)' 
Vnm/ 



(69) 



The near field condition requires that the distance between the atoms be much smaller than the wavelength associated 
with their first optical resonance. For hydrogen this wavelength is A = c/fi = 122nm. So, for the case where the 



prefactor of ( 69 ) is order unity and the interatomic distances are in the range of a few nanometers we find the 
entanglement force dominates over the standard London form see Fig [2j. 



Conclusion 



In this paper we have laid down the theoretical groundwork for the study of interatomic forces under fully noncqui- 
librium conditions. As a first step we have employed the influence functional formalism to derive a fully dynamical 
description of the atom-atom force for general atomic motion and initial states. We have found that a careful treatment 
of the infinite time limit shows the existence of a novel far field scaling when the atoms and field are not in thermal 
equilibrium. The dominance of this term in the laboratory would require a careful balance between temperature 
and the first optical resonance of the atomic species. For entangled atoms a novel near-field scaling is obtained that 
dominates the standard London force in certain regimes. These new forces could play an important role in quantum 
computing schemes involving entangled atoms. 
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